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Abstract 

We study the mean-field thermodynamic limit for a class of isolated 
Newtonian N-body systems whose Hamiltonian admits several addi- 
tional integrals of motion. Examples are systems which are isomorphic 
to plasma models consisting of one specie of charged particles mov- 
ing in a neutralizing uniform background charge. We find that in the 
limit of infinitely many particles the stationary ensemble measures with 
prescribed values of the integrals of motion are supported on the set 
of maximum entropy solutions of a (time-independent) nonlinear fixed 
point equation of mean-field type. Each maximum entropy solution of 
this fixed point equation can in turn be either a static or a stationary 
solution for the entropy-conserving Vlasov evolution, or even belong 
to a one-dimensional orbit of maximum entropy solutions which evolve 
into one another by the Vlasov dynamics. In short, the macrostates 
of individual members of an equilibrium ensemble are not necessarily 
themselves in a state of global statistical equilibrium in the strict sense. 
Yet they are always locally in thermodynamic equilibrium, and always 
global maximizers of the pertinent maximum entropy principle. 
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1 Introduction 



A short while ago I had an epiphany regarding equilibrium statistical mechan- 
ics. It happened while collaborating with Carlo Lancellotti on some curious 
Hamiltonian systems with unusual many body dynamics discovered by the 
Lynden-Bells |LBLB99| ILBLB04] . It could (and should!) have happened much 
earlier, though it wouldn't have been an epiphany, then, if I had understood 
right from the beginning what equilibrium statistical mechanics is all about. 
In any case, until that moment of revelation I was under the erroneous - 
yet I believe it's fair to say: common — impression that equilibrium statisti- 
cal mechanics is entirely about the statistical microscopic foundations of the 
macroscopic concept of thermal equilibrium for Hamiltonian iV-body systems 
when N 3> 1. Of course, equilibrium statistical mechanics is also about this, 
and even was invented specifically for this purpose, but that's not the whole 
story. The complete one is much richer and much more interesting. 

In the following, I will first recall the "folklore" about Boltzmann's |Bol96] 
ergodic ensemble measure, more familiar under the name Gibbs' |Gib02j mi- 
crocanonical ensemble measure 

d ^ = i^^(#-£) d3 V 3 \ (i) 

where 

$W(E) = j 8{H- E)d 3N pd 3N q (2) 

is the normalizing factor. I will point out misconceptions in this folklore which 
become particularly troublesome when considering the restrictions of ([1]) to its 
stationary ergodic sub-ensembles associated with additional conserved quanti- 
ties beside the Hamiltonian H. The true significance of ([1]) and its stationary 
sub-ensembles will become crystal clear once we pay due attention not just to 
the measures but also to the flow with respect to which they are stationary. 

Throughout this presentation, with the exception of a few more general 
results it is always understood that the Hamiltonian H of the classical N 
body system under consideration is stable [i.e., collapse of parts of the system 
to a point is excluded] and confining [i.e., escape of particles from the system 
is excluded]. Prime examples of such systems with long range interactions 
are one-component Coulomb systems in a neutralizing background, which in 
particular includes the Lynden-Bells' system (a 4D Coulomb system restricted 
to 3D motions) but also 2D point vortex systems. Yet, for simplicity here I will 
only present a general discussion for N body systems in 3D. The applications to 
the Coulomb systems will be published jointly with Carlo Lancellotti [KiLa08j . 

Finally, to simplify notation I absorb normalizing factors such as ([2]) into 
the measures by writing £(•••) for normalized delta measures, and I write the 
sets on which the 5's are supported as subscripts. 
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2 Folklore 



Without attempting to give a verbatim quotation from any particular book, 
I believe the following statement (S) quite accurately captures the essence of 
similar statements in the pertinent literature: 

(S) u An isolated classical N-body system with Hamiltonian H and energy E 
will spend most of its time in thermal equilibrium. In statistical mechanics, this 
equilibrium state is given by Boltzmann's stationary ergodic ensemble measure 



unless additional constants of the motion exist; in this case there is a stationary 
probability measure for each ergodic component M, given by 



where M = n£L {lfc = C^}, the Ik are the isolating integrals, with I = H, the 
Ck their values, and d 3N pd 3N q stands for Liouville measure on phase space." 

Since evaluating in one of the usual ways (that is approximately in 
most cases) will yield very satisfactory agreement with the empirical data al- 
most anytime we inquire into the thermal equilibrium properties of a many 
body system which can be described classically, and since the general em- 
pirical relations of thermal equilibrium are recovered rigorously from Q in 
some appropriate thermodynamic limit N — > oo dictated by the Hamiltonian, 
most of us, presumably, do not take any offense at reading statements such as 
(S), although we should. Indeed, in the line above ([3]), (S) makes the unfortu- 
nate suggestion that the macroscopic thermal equilibrium state of an individual 
physical system should be identified with a stationary probability measure over 
an ensemble of such systems, and this leads to some grossly incorrect conclu- 
sions in all but the most favorable circumstances. 

What's more, (S) is conceptually absurd. For instance, many invoke infor- 
mation theory to justify probability measures such as (j3J) and (@J, becoming 
subjective, and then (S) violates what should be enshrined as one of the 10 
commandments of physics: 



But even if the ensemble probabilities are interpreted in an objective man- 
ner the identification of the equilibrium macrostate of any individual physical 
system with a stationary ensemble probability measure is conceptually absurd, 
although in favorable cases when the overwhelming majority of systems in the 

1 Read: Do not confuse what a system does with what you can know about it! 
(The negation of this commandment is being worshiped in the orthodox quantum commu- 
nity, which seems exclusively concerned with observation and, most recently, information. 
But I don't want to take off on this tangent here. I have expressed my views in [Kie08| .) 



dfi = 5r H=E} (d m pd m q), 



(3) 



d/i = 5 M (d 3 V 3A g), 



(4) 
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ensemble have essentially the same macrostate then the results computed with 
this macrostate agree quantitatively with the results computed with the en- 
semble measure. Clearly, in such favorable cases some law of large numbers 
holds for the ensemble, i.e. if one formally let's N —>■ oo the empirical statis- 
tics of almost all members of the ensemble agrees with its ensemble mean. 
Boltzmann showed that such a law of large numbers holds for the perfect gas, 
assuming that H is the only conserved quantity and using the ensemble mea- 
sure ([3]). This prototype result by Boltzmann may have suggested the kind 
of thinking expressed in (S); all the same, (S) is absurd, and while masked 
by quantitative insignificance whenever a law of large numbers holds for the 
ensemble, the absurdity is there to stay and becomes manifest whenever the 
ensemble as a whole does not feature a law of large numbers — which is the 
case, for instance, at phase transitions of first order. Such phase transitions 
are defined quite sharply in the formalism of thermodynamics but not for ([3]) 
when N < oo, yet the situation is quite easy to understand: if roughly fifty 
percent of the ensemble members have "spin +1," the other half have "spin 
— 1" (to use a simple representative example to illustrate the idea), the ensem- 
ble mean is "spin zero;" yet, each member of the ensemble has either spin +1 
or spin —1, not spin « 0. 

Incidentally, Boltzmann was well aware of the fact that the notion of a 
phase transition and other thermodynamic notions as well are not unambigu- 
ously defined in statistical mechanics for any N < oo, and he also understood 
that all these thermodynamic notions can be sharply recovered from statistical 
mechanics in some limit N — > oo. To rigorously implement this idea took a lit- 
tle longer, until Ruelle, Fisher, Lebowitz, Lieb, Thirring and others published 
their seminal papers on the subject in the 60s and 70s of the previous century, 
see |Rue69l IThi89] . The insight that the thermodynamic notions are not exact 
for N < oo yet rigorously recovered when N — > oo has lead to the suggestion 
that statements like (S) have to be understood as implicitly referring to the 
(or some) "thermodynamic limit" N = oo, with a system's macrostate being 
identified in general not with the limit of the measure (djj or ^ itself but with 
the extremal ergodic measures into which the limiting ensemble measure de- 
compose^ — the so-called pure states. This avoids the spin-up vs. spin-down 
problem of the above example (and all problems with similar dichotomies): 
the extremal ergodic measures are precisely the spin— 1 and spin+1 measures. 
The pure states of the N = oo ensemble define ergodic sub-ensembles; a law of 
large numbers holds for these sub-ensembles though not for the full ensemble. 
But has this now eliminated the conceptual problems? 

Although the invocation of the limit iV — > oo combined with the decom- 
position into extremal measures avoids the problems with dichotomies, the 

2 The set of probability measures on some measure space is convex, and every point in a 
convex set can be written as a convex linear combination of the extremal points of the set. 
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proposal to define what is meant by a (macro)state of an individual iV-body 
system (as distinct from computing it to a high degree of precision) in terms of 
the limit N — > oo was once commented on by Wigner thusi " [Name], you are 
a great man, but you are are not infinite." Wigner's point, translated into our 
setting, was that even granted the extremal ergodic measures of the N = oo 
ensemble were legitimate macrostates of the N = oo system (they are, as we 
shall see, though there is still a subtlety to be addressed), that still doesn't 
free us from the obligation to define a macrostate of a finite- iV system without 
invoking the limit N — » oo. Of course, having found a sensible definition of 
macrostate for the N = oo system, a finite-iV definition should coincide with 
the N = oo definition in the limit. I will give a legitimate definition below. 

I now address the second distinction of importance to which I have alerted 
you by highlighting in the paragraph after (S): the distinction between the 
equilibrium of an individual system and the stationarity of an ensemble mea- 
sure. It is very easy to be confused about this distinction, even if yon are not 
confused about the distinction between macrostate of an individual system 
and the state (i.e. probability measure) of an ensemble of such systems as 
discussed above. Part of the confusion is just a matter of semantics. With 
hindsight, Gibbs' choice of terminology, to call ensembles defined by ([3]) or 
(|4j) statistical equilibrium ensembles |Gib02] . while sensible, was not very for- 
tunate. Another semantic item which contributes to the confusion is the fact 
that thermodynamic equilibrium of an individual iV-body system has become 
synonymous with its statistical equilibrium, which terminology is based on the 
same notions which are behind Gibbs' terminology for a stationary ensemble 
of systems, only now applied to an individual system viewed as an ensemble 
of particles. The provocative title of my presentation is deliberately chosen 
to put you up against all the confusion right at the start: I insist, the title 
makes sense — but which sense? To find out we have now to confront the 
serious confusion, which is not at all a matter of semantics but belongs in a 
category described beautifully by Franz Kafka, who said "To correctly grasp 
an issue and to misunderstand the same [issue] do not exclude one another." 
(Translated from a quotation in |Rue69j .) 

A stationary ensemble measure by its very definition does not depend on 
time, and this of course remains so also when taking the limit N — > oo. The 
limiting measure has a unique convex decomposition into extremal measures, 
and since the limiting measure is time-independent, so is its extremal decom- 

3 As told to mc by Shelly Goldstein. Wigner commented on the suggestion by a distin- 
guished seminar speaker that the measurement-induced collapse of the n-body wave function 
of the measured system, postulated in orthodox quantum mechanics, would perhaps result 
rigorously from the unitary Schrodinger evolution of the combined, n + TV-body system in 
the limit N — > oo, just as sharp phase transitions do. Clearly the gist of Wigner's remark 
is also valid in our context. 
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position. If these extremal states represent the N — > oo limits of the possible 
macrostates of individual finite- iV systems in the ensemble, as we said earlier 
they would, only time-independent macrostates can be obtained. 

Here then is something paradoxical: all the measures (jlj) are stationary, 
so what I just said applies to them. In particular, it applies to the special 
ensemble measure of the type (j3J) which is the next best known after ([3]), 
namely the one for M = {H = E} fl {J = L} where the Hamiltonian H 
and the z-component (say) J of the angular momentum are conserved, having 
values E and L, respectively. Now when L ^ the system must rotate, in 
fact rigidly as Landau and Lifshitz will tell you. But we know that there are 
situations where the rotating system is not rotationally invariant, and then the 
macrostate of this system will be time-dependent, of the kind cos(cot— ip), with 
uj the constant angular frequency of the rigid rotation. For example, think of 
the triaxial Jacobi ellipsoids which describe rigidly rotating constant density 
self-gravitating fluids |Cha69] . How does this all fit together? 

According to the folklore the paradox is resolved as follows: 
(S') " Since a stationary ensemble only yields stationary macrostates (read: 
equilibrium states), and since the rigidly rotating configurations with non- 
vanishing angular momentum rotate with a constant angular frequency, one 
has to apply the methods of equilibrium statistical mechanics in a co-rotating 
frame where the macrostate of the individual system is now stationary." 
With regard to an inertial frame such rotating macrostates are usually called 
rotating equilibria. This prescription of how to apply the equilibrium statistical 
mechanics formalism leaves a strange taste in one's mouth, for it avoids the 
answer of the puzzle by circumventing it; yet the described procedure has 
worked fine since its inception, also for me. 

After relegating the above paradox to the back of my mind for many 
years, it finally came full back into my conscience after Carlo Lancellotti and 
I had begun to work on the rigorous derivation of what the Lynden-Bells 
|LBLB99| ILBLB04] had called a "perpetually pulsating equilibrium." They 
had discovered a curious Hamiltonian iV-body system in three dimensional 
space with an unusual additional conservation law, besides the conventional 
ones for energy, momentum, and angular momentum. They managed to show, 
by studying the virial theorem, that the mean-square radial distance from 
the center of mass undergoes a sinusoidal motion for all N. They had also 
carried out some computer simulations with moderate N and looked at the 
empirical statistics of velocities in a "co-pulsating frame," finding evidence for 
a Maxwellian — hence their terminology. But when Carlo and I set out to 
apply the conventional strategy advocated in the folklore (S'), i.e. transform- 
ing the problem into a co-moving frame, I soon grew dissatisfied with this 
approach for conceptual reasons. The dynamics of the co-moving frame, while 
still simple, is already quite non-trivial because now the angular frequency 
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of the rigid rotation is itself time-dependent (think of the familiar sight of 
an ice skater performing a pirouette, with her arms periodically extending 
from and retracting to her body). If we permit terms like "rotating equilib- 
rium" and "pulsating equilibrium," we certainly have to also allow "rotating- 
pulsating equilibrium," as this is what the Lynden-Bells' "equilibrium" does 
in general. But what if someone discovers an even more curious Hamiltonian 
iV-body system, with a "rotating-pulsating-twisting-bending-shearing-jerking- 
and-what-have-you equilibrium" — where does it stop? Are we not diluting 
the meaning of the word "equilibrium" until the construction becomes totally 
absurd? 

Somehow the folklore was totally off target. 

3 Epiphany 

The crucial insight is that something everyone knows to be true anyhow for 
the finite- N measures (jlj) in regard to the actual finite- N microstates is true 
- by inheritance — also for the macrostates of the actual individual systems 
in the ensemble, and remains true in the limit iV — > oo. 

In a nutshell, (jlj) lives on the iV-body phase space, which consists of the 
generic phase points of iV-body systems. Any point in iV-body phase space 
defines a possible microstate of an iV-body system. In addition, there is a 
Hamiltonian flow on this generic iV-point phase space. Generic points do not 
move, they are independent variables. What evolves in time is the actual point 
in phase space which identifies the microstate of the actual iV-body system. Its 
motion follows a particular flow line (if the Hamiltonian is time-independent) 
picked out by the initial actual phase point. So while (jlj) is stationary (with 
respect to the flow on iV-body phase space), the ensemble of actual iV-body 
systems it represents consists of systems whose actual phase points in general 
will not be stationary. Everyone knows that, and no one has a problem with 
that, even though some people may prefer to say the same thing differently. 

Now, any mapping of the iV-body microstates into iV-body macrostates 
also induces a mapping of the flow on the space of generic microstates into 
a corresponding "flow" on the space of generic macrostates. Let's call this 
new "flow" the macroflow. Also the measure (jlj) will accordingly be mapped 
into a stationary measure on the space of generic macrostates. Let's call that 
induced ensemble measure the macrostate measure. If in the limit N — > oo the 
macrostate measure is supported only on such generic macrostates which are 
stagnation points for the limiting macroflow, then the limiting macrostates of 
the actual systems will be stationary too, and so the true finite- N macrostates 
will be approximately stationary, with small fluctuations. It is in the nature 
of macrostates that they don't care about most of what is going on at the 
microscopic level, so macrostates can be relatively quiet even if the underlying 
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microstates evolve violently. But if that limiting measure is supported on 
generic macrostates at which the limiting macroflow is not stagnant, then 
the limiting macrostates of the actual systems will be dynamical, traversing 
an orbit in the subset of generic macrostates which form the support of the 
limiting stationary macrostate measure. By the same token as before, the true 
finite- N macrostates will exhibit similar large scale macroscopic motions. 

It should be perfectly clear by now that there is no conflict between the 
ensembles being stationary, implying that their extremal decompositions in 
the limit N — > oo are stationary as well, and the undeniable fact that the 
macrostates of individual actual systems making up the ensemble will in gen- 
eral exhibit macroscopic evolution in time. 

In the remaining sections I will carefully explain how this general story pans 
out for the ensembles (j3J). I hope that my explanations are entirely correct, 
but I cannot rule out that I am still confused about this or that — think of 
Kafka! As appropriate for this conference, the general story is worked out for 
systems with long-range Galilei-invariant pair interaction which I will specify 
in more detail below; yet, the next section and part of section 5 are valid for 
Hamiltonians H with more general Galilei invariant pair interactions. I will 
not invoke transformations into any co-moving non-inertial frames whatever; 
however, a Galilei boost into an inertial frame in which the center of mass is 
at rest is of a different category and will be performed (though the formalism 
does not even require that!). Of course, some limit iV — > oo will be invoked, 
though not to define macrostates of the finite- N systems but to calculate them 
approximately. Macrostates will be defined for finite- iV systems without invok- 
ing that limit. In fact, our first duty is to define what we mean by microscopic 
and macroscopic states, or micro- and macrostates for short. There is also an 
intermediate level frequently referred to as mesoscopic. 

4 Micro-, meso-, and macroscopic 

Taken literally the term "...-scopic" means something like "point of view," a 
very subjective notion, and we sure don't want to confuse what systems do 
with how we look at them! Fortunately there is a different way of putting 
it, namely as the objective structures produced by the iV-body dynamics at 
the "finest scales," at "intermediate scales," and at "large scales." There is 
still some ambiguity as to what exactly is "fine," "intermediate," and "large," 
and this ambiguity will not go away! Yet at least the reference to subjective 
notions which seem implied by "scopic" is gone. 

However, concepts must not be too arbitrary or they lose their significance. 
Since by "state" of a system one means a complete specification of the dynam- 
ical variables of the system, it must evolve (to good approximation) according 
to an autonomous dynamics at the scale under consideration] cf. |GoLe04j . 
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4.1 Microstates of individual systems 

Although the term microscopic scale does contain the ambiguity aluded to 
above, there fortunately is a clear cut choice for what the microscopic state of 
a Hamiltonian iV-body system should be. Our Hamiltonian systems consist of 
N > 2 interacting nonrelativistic Newtonian point particles in Euclidean E 3 
with generic momenta p i G R 3 and positions q i G M. 3 in an arbitrarily chosen 
inertial frame 1R 3 . Thus, for our Hamiltonian iV-body systems then any point 
(p 1 , ...,p N , q 1: q N ) = X G M 6Ar represents a possible state of the iV-body 
system, or a generic state. If q^t) and Pi(t), i = 1, N, are the N particles' 
actual position and momentum vectors at time t G M in an inertial frame, then 
(Pi(t), ...,q N (t)) = X t is the actual state at time t. This is the most detailed 
characterization possible for the Hamiltonian iV-body system, and therefore 
this is its microstate. This is the conventional definition. It is natural and 
unambiguous. 

The disadvantage of this conventional definition of the microstate is that 
it is not well tailored to our needs. So instead we will work with a completely 
equivalent characterization of the microstate of a Hamiltonian iV-body system, 
equivalent in the sense that there will be a one-to-one mapping between both 
types of objects. We identify each point (p 1 , q N ) = X with a singular 
empirical "density" (a measure) on M 6 as follows: 



Obviously the map is bijective, i.e. X <-> Ax- Henceforth we will refer to ([S]) 
as the generic microstate of an iV-body system. Note that ([5]) means physically 
the number density of particles at q G R 3 with momentum p G M 3 . Since by X t 
we denote the actual point in phase space R 6iV at time t, then for us Ax t (p, q) 
is the actual microstate of the iV-body system. 

For the sake of completeness, we remark that this concept of microstate 
extends also to an oo many particles system. Intuitively: countably infinitely 
many points in M 6 have Lebesgue measure 0. Even though the empirical 
measure YHfLi $(P ~ Pi)^{Q ~ Qi) * s n °t even locally finite if oo many points 
are densely distributed in some domain C M 6 , one can make sense out of it by 
looking at it as the collection of its finite iV truncations equipped with total 
variation (TV) topology. In TV topology then there will be no convergence 
to a continuum function when iV — > oo, so even though the oo many points 
may be dense in some domain, the topology is fine enough to distinguish the 
state from a continuum state. However, we will not need this construction. 



N 




(5) 
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4.2 Mesostates of individual systems 

We now face an irreducible ambiguity: what exactly constitutes a "mesoscopic 
scale"? There is no such thing as "the mesostate" of an iV-body system per se. 
Yet in statistical physics one usually invokes this terminology in the context of 
kinetic theory, when a system is described in terms of some continuum formal- 
ism in (p, q) space M 6 ; think of Boltzmann's transport equation for a dilute 
gas, or Vlasov's equations. Accordingly we could be tempted to call any such 
continuum density function f(p, q; t) of kinetic theory an "actual mesostate," 
and the time- independent generic functions the "generic mesostates." This is 
not quite adequate, though. The reason is that a continuum theory is of course 
even more detailed than an iV-body theory, since / needs to be specified at 
uncountably many points. This is hardly a description that considers only 
scales sufficiently bigger than the microscopic one of the previous subsection. 

This inadequacy of the continuum description is ameliorated by restrict- 
ing the (positive) continuum density functions f(p,q) in a suitable way. For 
instance one could work with piecewise linear, Lipschitz continuous functions 
with Lipschitz constant < C for some C, and which linear pieces have their 
epicenters at a fixed "hypercubic" lattice. This would be quite adequate. 

To make our life easier we are a little more cavalier and allow functions 
f(p,q) which are uniformly Lipschitz continuous and uniformly iterated Lip- 
schitz continuous. Basically this says that such functions cannot be too steep, 
and cannot oscillate on too small scales; and since they are non-negative and 
need to integrate to a finite number, they cannot develop singularities when 
(p, q) go to infinity. We denote this class of functions j^S '. It is understood 
that some uniform bound on the Lipschitz and iterated Lipschitz constants is 
chosen; they clearly restrict the possibilities, but how exactly is a matter of 
choice and good judgement and part of the irreducible ambiguity of any notion 
such as mesoscopic, or for that matter also of macroscopic. 

Next we associate to each microstate Ax a mesostate fx- This map may 
be but need not be bijective. To simplify the discussion, we work with the 
normalized empirical densities A x = N _1 Ax and understand henceforth that 
any fx is normalized such as to integrate to 1 as well. We now pick a suitable 
Kantorovich-Rubinstein metric oIkr (which metrices weak convergence of nor- 
malized measures) and define fx as the minimizer of the KR distance between 
A x and i.e. 

^kr (A x , fx) = mm{d KR (4 X , /) |/ e JZ£} (6) 

if the minimum exists; if it doesn't, since the infimum of the set on the r.h.s. is 
certainly strictly positive, one can stipulate some convenient rule to choose an 
fx among those / which differ from A x in KR distance by not more than the 
infimum plus e. If that sounds somewhat ambiguous, that's because it is, but 
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such ambiguities are acceptable given the inherent ambiguity in the notion of 
a "mesoscopic scale." 



4.3 Macrostates of individual systems 

Of course, the same kind of irreducible ambiguity that haunts the "meso- 
scopic scale" also raises the question of what exactly constitutes a "macro- 
scopic scale;" and again there is no objective answer to this. In the statisti- 
cal physics literature, the main distinction between the "macroscale" and the 
"mesoscale" is that in the former one now does not resolve the variations in 
p space. More to the point, a macrostate is usually the small collection of 
fluid dynamical functions on q space M 3 , like mass density, velocity field, and 
kinetic energy density field or pressure field. Such functions can be obtained 
from the mesostates by taking the first few moments over p space. 



5.1 Dynamics 

The dynamics of each system is governed by a Hamiltonian of the type 



H [N \p 1 ,q 1 ]-..]p N ,q N ) := ^ 2^7 \Pi\ 2 + Yl e ^ V (\ 9i ~ Qj |)> ( 7 ) 



where the > are the inertial masses and the e« G M "generalized charges." 
The pair interaction potential V(r) will be of long range and locally integrable 
to admit a mean- field limit when iV — > oo, but many facts do not depend 
on all the details of V. In particular, the following familiar conservation laws 
hold for a larger class of V than just those allowing mean-field limits. Since 
we need the conservation laws to set up our ensembles, we collect them here. 

Lemma 1. Assume the pair interaction potential V satisfies V G C 2 (1R + ,M) ; 
guaranteeing local existence and uniqueness of the dynamics. Then by the in- 
variance of under time and space translations, as well as space rotations, 
the isolating integrals of motion are: the Hamiltonian given in Ffy, the 

total momentum phase space function 



5 Finite N 



Ki<N 



l<i<j<N 





Ki<N 



and the total angular momentum phase space function 



(9) 



Ki<N 
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Lemma 1 exhausts the isolating integrals of motion w.r.t. the center-of- 
mass frames which are shared by all Hamiltonian systems characterized by (171) 
under the stated hypotheses on V. Additional isolating integrals for all N > 2 
may exist for special cases of V, and we shall encounter an example later on. 

Moreover, on account of H^ N ' changing only by an additive constant un- 
der a Galilei boost, what is sometimes called the system's centroid and given 
by ]^Ei<i<iv( m ^i(*) ~tPi(t)), is conserved, too; here M = Ei< t <v m * is 
the total mass of the system, which is independent of the choice of inertial 
frame. However, the centroid is not the evaluation of a t-independent func- 
tion on (pi, . . . , q N ) space and therefore not an isolating integral. Yet, since 
by the Galilean invariance of Newtonian point mechanics any particular one 
of our iV-body systems can always be described w.r.t. any of its center-of- 
mass frames (any co-moving inertial frame for which the center of mass is the 
origin of that frame), without loss of generality it suffices to study iV-body 
systems under the holonomic scleronomous constraints P ( - N \p l , . . . , q N ) = 
and C {N \p l ,q 1 ;...;p N ,q N ) =0, where 

C (N) (p 1 ,q 1 ;...;p N ,q N ) := — ^ m iq% . (10) 

l<i<N 

is the center- of -mass phase space function. As far as the stationary ergodic 
subensembles go that we introduce below, these center-of-mass-frame con- 
straints can be treated on the same footing as the the true isolating integrals. 

Hence, to unify the notation, we include both the true isolating integrals 
and the scleronomic constraint functions in a list of functions I k , k = 0, 1, 
simply called invariants, with C k denoting the values they take. In particular, 
we define: I := # {A ° with C = E, (h,h,h) ■= P {N) with (Ci,C a ,C 3 ) := 
0, (h, h, h) := J {N) with (C 4 , C B , C 6 ) = L, and (J 7 , J 8 , J 9 ) := C (7V) with 
(C*7, Cs, Cg) = 0, where E G K and Lei 3 are time-independent parameters. 

Incidentally, the phase space functions ©, ([8]), Q), (|T0|) are not all in 
involution; the invariants H^ N \ P^ N \ {J^l 2 , ■ e are, where e is any 

fixed unit vector in space. Moreover, in any of our center-of-mass frames 
is in involution with but not with p( N > and J^ N \ This is not a problem, 
for we will characterise the ergodic submanifolds conveniently in such a way 
that we do not need to find a set of invariants in involution; however, it is 
implicitly understood in the following that any additional invariants Ik with 
k = 10, ... will not be redundant. 

Lemma 2. Let K + 1(> 10) be the number of invariants admitted by H^ N \ 
with V as in Lemma 1. 

a) Let WIk denote the manifold 

M K = {( Pl ,...,q N )e R m \I k = C k , k = 0, K) . (11) 
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Then M.k is invariant under the flow generated by the Hamiltonian . 
b) Let fj,^ denote the singular measure supported on M# given by 



(12) 



Then [W\l is invariant under the adjoint flow associated with H^ N \ i.e. [W\l 
is a stationary weak solution of Liouville 's equation 



So far, the conditions for Lemmas 1 and 2 also cover Hamiltonian systems 
of type (JTj) with physically relevant pair interactions, such as V(x) = K/x with 
Ci G M. (Coulombian electricity) and V(x) = —K/x with ej G M + (Newtonian 
gravity), as well as V(x) = K(x~ 12 — x~ 6 ) with G M. + (Lennard- Jones 
molecular pair potential). Unfortunately, neither of these Hamiltonian systems 
is self- confining, and the jury is still out on the question of the generic global 
existence of the dynamics in the gravitational and electrical iV-body problems. 
We now impose stronger conditions on V which guarantee global existence of 
the dynamics and self-confinement of the system of particles. 

Proposition 1. Let K + 1(> 10) be the number of invariants admitted by 
H( N \ with V satisfying the hypotheses in Lemma 1. 

a) If in fact V G £ 2 (M + ,M+) with -e^V^r) ^ oo for all or V G 

£ 2 (IR + , M + ) with V'{r) ; and in either case lim sup^^ |V(r)|/r < C, 
then the dynamics exists globally in time. 

b) If in addition to the hypotheses in a) we have lim^oo eiCjV{r) = +oo 
for all pairs i,j, then [W\l is a finite measure. 

Henceforth we assume that V satisfies the hypotheses of Proposition lb). 
Note that when N > 2 this implies that either > Vz or < Vi. 

5.2 Statistics 

In the section "Folklore" I spoke several times of "a law of large numbers" 
as iV — > oo. This is related to, but not to be confused with, a familiar law 
of large numbers which holds (under mild conditions) for ensembles of i.i.d. 
finite-iV(!) systems. 

Denoting normalized (fl2l by i.e. f d//^ = 1, this can be interpreted 
as the single-system a-priori probability measure for a Newtonian iV-body 
system in phase space M. 6N to be at X if all that is given about the system is 
its Hamiltonian H^ N ' with energy E, its total angular momentum L, its center 
of mass = 0, its total momentum = 0, and the values Ck for whatever other 
invariants there are; with also the number iV of particles and total mass 



d tf i+{ f i,H^} = 0. 



(13) 
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M are given. Thus, if S3 denotes the Borel sets of R 6N , then (R 67V , S$^ N) ) is 
our single-system probability space. The conventional microstate of a single 
system is now a random vector X = (CPi , . . . , JV, Qi, . . . , Qn) G K 3Ar x R 3N , the 
specification of which amounts to giving exactly the positions and momenta 
of each particle in the iV-body system. If B C M 6N is a Borel set, then the 
probability of X being in B is 

Prob (X G B) =^ N \B). (14) 

Clearly, Prob (X G B) = unless B HM K ^ 0. 

The Boltzmann ergodic subensemble associated with f|T2|) is an infinite 
family of i.i.d. random vectors X^, each with a-priori distribution (fT2l . and 
each one of which represents a Newtonian iV-body system which is governed 
by the same Hamiltonian H( N ' and has the same energy E, center of mass 
= 0, total momentum = 0, total angular momentum L, and same values Cf. 
for whatever other invariants there are; of course, each system also has the 
same number N of particles and total mass M. Any ordered set of £ copies 
of X is a random vector (X^\ ...,X^) G M. i6N with joint distribution ju^'^. 
For this family of i.i.d. systems the conventional weak law of large numbers 
(WLLN) holds. 

We are going to state the WLLN in the form adapted to our discussion 
of the concept of a system's "microstate" as A x (p,q)- Accordingly, X is 
identified with a singular-empirical-density- valued random variable Ax(f», q), 
and the measure (1T21) is identified with a single-system a-priori probability 
measure jl^ on the convex set of all singular densities (i.e. measures) on IR 6 . 
Clearly, /}W is supported only on the set of all singular densities of the type 
Ax(p,q)- Boltzmann's ergodic ensemble is therefore identical to an infinite 
family of i.i.d. random densities Ax(p, q), each with a-priori distribution p,^'. 
By convexity, the empirical mean of a size-£ sample of i.i.d. empirical densities, 

A ( x ) ( P ,q) = ~/£A xU) (p,q) (15) 

3=1 

is again a density of the same total measure = N. For the normalized measure 
fj,( N > defined by (fi~2l) . denote by /i^p (d 3 pd 3 g) its first marginal measure for the 
z-th particle variables (p i( qrj, evaluated at (p,q). Recall (see, e.g. [Dur96j 
IDud02j ) that a family of probability measures v n on M. d is said to converge 
weakly to v if f g(x)dis n — > f g(x)du for all bounded continuous functions 
g G £°(IR d ). Recall furthermore (e.g. [Dud02j) that weak convergence is 
metrized by some Kantorovich-Rubinstein metric g?kr- We are now ready to 
state the weak law of large numbers. 
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Theorem 1. For Boltzmann's ergodic ensemble of i.i.d. random vectors X, 
each with a-priori distribution fHfy. we have the weak law of large numbers, 

Jim Prob (d KR (A?(p l9 )dV 3 g, E^M?(d 3 pd 3 ?)) > c) = 0. (16) 
We follow the probabilists custom and rephrase the WLLN shorter thus: 

JV 

A%(p, q)d 3 pd 3 q £ ^ /x 1 .(d 3 pd 3 g) in probability. (17) 

i=i 

So far all the particles in an individual iV-body system may have different 
masses and chargesQ Many systems of interest in nature do consist of only a 
handful of species, and the prototype systems to study are one-specie systems 
with rrii = m > and e, = e G K for all % in ([7]), so the Hamiltonian becomes 

H^( Pl ,..., qN )= y: iw 2 +EE eV (k-^D- ( 18 ) 

l<i<N l<i<j<N 

For such single specie systems the WLLN simplifies to the following. 



Theorem 2. For the Hamiltonian < f7oj) . all the first marginal measures are 
identical due to the permutation invariance of which is inherited by 

the product over all the isolating invariants in fU| ). Thus, fi ( ^p(d 3 pd 3 q) = 

/j,t N \d 3 pd 3 q) for all i = 1, N. Accordingly, as £ — > oo, 

A^ip, q)d 3 pd 3 q — iV^)(dV 3 g) (19) 
m probability. 

Note that no limit iV — ► oo is involved at this level. 



6 The limit N — > oo 

From now on our Hamiltonian will be (TT5T) . In addition, we are interested now 
only in those Hamiltonians (Tl8l) which permit a mean-field continuum limit 
iV — > oo. Therefore, in addition to the hypotheses of Proposition lb, from 
now on we assume also that V satisfies V o | . | e ^^(IR 3 ), and we also ask 
that lim^oo V{r)/r = oo. 

4 Physical examples of large ./V-body systems with a distribution of particle masses are 
stellar clusters. 
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6.1 Dynamics of microstates 

Having identified points X on C M 67V with empirical densities Ax(p, q), 
we note that the Hamiltonian flow on M# which transports any point Xq 
on M.k at time to a corresponding point X t on at time t also defines 
the dynamics of the empirical densities Ax{p, q)', namely, the initial empirical 
density A Xo (p, q) evolves into A Xt (p, q) at time t. Under suitable assumptions 
the dynamics of the normalized empirical densities A Xt (p, q) converges with 
N — > oo to a limiting dynamics of normalized continuum densities f(p, q; t). 

In the following, we attach a superscript ( N > to M K and X 6 to indicate 
that M K = M^ and X = X™ change with N. We also write ijW(I^) 
for H^{p 1} ...,q N ), and similarly for the invariants P^ N \ C^ N \ and J^ N \ 

We are interested in a continuum limit in the sense that the particle po- 
sitions q { (we should write □;) converge in distribution as N — > oo, i.e. (at 

time t = 0) we want / R3 A_ x(N) (p, . )d 3 p p( . ; 0) G n C 6 °)(M 3 ) weakly 
(metriced by some Kantorovich- Rubinstein metric c?kr)- We are only inter- 
ested in those p for which JfV(\q — q'\)p(q;t)p(q';t)d 3 qd 3 q' < oo, where JJ 
is taken over M 3 x M. 3 . We also want the kinetic and potential energies to 
scale comparably with N for large N. Since only N summands contribute 
to the kinetic energy while the potential energy oc N 2 (for N ^> 1) when 
J R 3 A X ( N ) (p, q)d 3 p — > p(q) in the stipulated sense (suppressing t temporar- 
ily), to achieve equal iV-scaling of kinetic and potential energies the rescaled 
particle momentum vectors N~ 1 ^ 2 p i must converge in distribution, so that 

N ' 2 Ei<i< N \Pi\ 2 2£ kin- Thus we ask that N 3 l 2 A x{N) (N l l 2 p,q) 

fiP, q) e (£+,i n C°)(R 6 ) weakly, such that N~ 2 H^(X^) N -^? e(/) = e < 
oo, where 

£(/) = ^JJ\p\ 2 f(p,q)d 3 pd 3 q 

+ j ffjf V(\q - q'\)f(p, q)f(p', q')d 3 pd 3 qd 3 p'd 3 q> , (20) 

is the "energy of /." It then follows also that| C iN) {X^) £(/) = 0, 

iV-3/2pW(xW) ^5° *(/) = 0, and iV" 3 / 2 J (7V) (X^)) j(f) = I, where 



oXf) = JJqf(p,q)<i 6 pd 6 q (21) 

*(/) = f[pf(P,q)d 3 pd 3 q (22) 



5 The vanishing of Q,(/) and 2>(/) follows of course from the fact that each is picked 

from a M^°, so that necessarily C W (XW) = and pWpfW) = for all N. 
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/(/) = |Jgxp/( Pl q)dV 3 g (23) 

are the "center of mass — ," "momentum — ," and "angular momentum of /," 
respectively. The functionals j(f), &(f), Q,(f) are well-defined for / G £,\ 1 
whenever £ (/) < oo thanks to our hypotheses on V. 

We are now ready to state our first dynamical mean-field limit theorem. 
It is understood throughout that V in (ITS]) satisfies the stipulated hypotheses. 
Recall that the list of invariants {Iq,...,Iq} is common to all Hamiltonian 
systems with a Hamiltonian of type (FTSl) . 

Theorem 3. Let Hg N ^ denote any Hamiltonian of type (T72J) for which the list 
of invariants {Jq, ...,/g} exhausts the invariants. Let 

% = {/ G 4 :1 (M 6 ) :£(/)< oo, |*(/)| < oo, < oo, |7(/)| < 00} (24) 

denote the set of admissible densities on M. e . Let {t 1— > -Xt^}jV=2 denote a 
sequence of solutions of the Hamiltonian dynamics on Mg , iV = 2,3, gen- 
erated by Hg . For the initial data, assume that N 3 ^ 2 A_ x (n) (iV 1//2 p, q) A ^5° 

f(p,q;0) G 2t 9 weakly. Then N 3 / 2 A (N) (N^p.q) /(p,g;t) G 2l 9 



-l/2 t 



weakly, where f(p, q; t) is a weak solution to the Vlasov equation 

d t f + m- l p • V/ - V (eV * / fd 3 p) -d p f = (25) 

with initial data f(p,q;0). Moreover, N' 2 H^(X^X/2 t ) -> £(/(.;*)) = e, 
- £(/(•;*)) = 0, and A^P^X^) - *(/(.;*)) = 0, 
and N~ 3 / 2 j( N \Xj**} 1/2 ) — ► /(/( • ; £)) = Z, wt/i £ one? Z independent of time. 

Beside the invariants //. inherited from the isolating invariants the 
Vlasov dynamics conserves an infinitude of so-called Casimir functionals of 
/■ 

Theorem 4. Let C : R + — > K 6e a continuous function for which c[f( ■ ; 0)] = 
// • ; 0))d 3 pd 3 g exists /or / G 2l 9 . Lett ^> f{.\ t) solve with initial 
data /(.;0). T/ien c[/(.;0)] = /or aZZt. 

In particular, among the Casimir functionals of (H3J] are: 
the "normalization of f, 

*(/) = jjf(p,q)d 3 pd 3 q, (26) 
£/je "entropy of f , " 

S{f) = ~II 9) ^ 9)d3pd3g • (27) 
6 Of course, #(/) = 1 for / G 2t 9 . 
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If hw given in (|T8l) is not of type Hg , then additional invariants 
(k = 10,. ..,K > 9) exist (because the list (I 1 ,...,I 9 ) is shared by all 
given by f[T8"j) . The Lynden-Bells' Hamiltonian is of type iijg ^ or which the 
obvious generalization of Theorem [3] holds, while Theorem H] continues to hold 
unchanged. 



6.2 The support of the stationary ensemble measures 

We begin by stating our mean-field limit theorem for // 7V - ) ; or rather: fx^ N \ 

Theorem 5. Let H {N) = H^ N) as in Theorem^ set E = N 2 e and L = N 3/2 l, 

with e > and I fixed, viable data. Let be given in [Tfy) with K = 9. 
Then 

lim = ^ (28) 

exists, and jx is a stationary measure for the adjoint dynamics inherited from 
H^ N ' in the limit N — > oo. Furthermore, jx is supported on the set of maxi- 
mizers of the entropy functional s (f) , given in (27\ ), taken from the set of trial 
densities 

X 9 = {/ G 4 5l (M 6 ) :£(/) = = 0, i(f) = j(f) = I}. (29) 

If f,(p,q) is a maximizer of s(f) over the set X 9; then f,(p,q) satisfies the 
fixed point (Euler-Lagrange) equation 

where p. = J f.d 3 p and where X^, A E; \ } , X e , A Q are the Lagrange multipliers 
for the constraints ot{f) = 1, £(/) = £, /(/) = ^ #(/) = 0, and Q,(f) = 0. 

We are primarily concerned with the support of fx. The interesting question 
of how fx varies over its support will be touched upon but briefly. 

In addition to the maximizers of s(f) there may be other solutions of (!30|) . 
yet only those solutions /, for which s(f») = max over Xg are in the support 
of jx. All solutions are critical points for s(f). Independently of whether a 
solution of (1301) maximizes s (/,) = max over Xg or not the following holds. 

Theorem 6. Every solution of the fixed point equation UJtyl satisfying the 
constraints ^C(f») = 1, &(f*) = ; &(/•) = ; and y(f.) = I factors into 
a product of a locally (at q) shifted Maxwellian on p space and a purely 
space- dependent Boltzmann factor with X T = 0, A E = — 1/k^T < 0, and 
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Xj = uj/k-QT, where T is the temperature and uj the angular frequency vector 
of rotation of the dynamical system; thus 

3 

q) = (^kr) 2 exp (-2^r lp ~ muJ x qf ) pM) (31) 

with p. satisfying the fixed point equation on q space 

exp (-jig, [e 2 (V*p.)(<?) - f |w x g| 2 ] + A a • q) 

P. (9) = k (32) 

/*> ex P ("5^ [e 2 (V* P .)(q') - f |w x gf] + A Q • g'j dg' 

Here we also eliminated the Lagrange multiplier using ^i{f,) = 1. 

Theorems [5] and [6] hold for all Hamiltonian systems with Hamiltonian of 
type Hg N \ which covers - presumably - the majority of models with Hamil- 
tonian given by (fl8|) satisfying the hypotheses of Proposition lb together with 
those stipulated at the beginning of the section on iV — > 00. They are proven 
similarly as in |KiLe 97j. It is straightforward to generalize Theorems [5] and 
[6] to the Hamiltonians given by ( fl8l) which are of type h[q , which includes 
the Lynden-Bells' Hamiltonian. In that case the entropy maximizer factors 
into a product of a locally (at q) shifted Maxwellian on p space and a purely 
space-dependent Boltzmann factor. The shifted Maxwellian in (13TT) is known 
as a "rotating Maxwellian;" in the case of the Lynden-Bells' Hamiltonian we 
find a "rotating-dilating Maxwellian." 



7 Thermodynamics [sic!] 

By Theorem [HI w.r.t. Kantorovich- Rubinstein topology (metrizing weak con- 
vergence of measures) every solution of the fixed point equation fl5U|) represents 
the microscopic density function on single particle phase space of an individual 
N = 00 many particles system which locally at each q is in thermal equilibrium 
w.r.t. an instantaneously co-rotating frame at q. A forteriori this is true also 
for the maximizers of s (f) over T9. However, globally a solution to ( 1301) — and 
therefore in particular any maximizer of s(f) over Xg — may be a stationary 
solution of fl25|) but need not be. If a solution to fl30|) is a stationary solution 
of ([25]) . then it describes a system which as a whole is either in a non-static 
stationary state or in a static state, depending on whether a non-vanishing ve- 
locity field u{q) = m" 1 J pf,(p, q)d 3 p exists or not, respectively. If a solution 
to ( |30l) is not a stationary solution of ( 1251) . then it is necessarily a snapshot of 
a whole connected orbit of densities which solve the genuinely time-dependent 
Vlasov equation ([25j) . Note though that all snapshot densities in such an orbit 
are necessarily maximizers of s (/) over Tg, for the Vlasov dynamics preserves 
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the isolating invariants as well as normalization and entropy. All this will be- 
come more dramatically so for the rotating-dilating Maxwellians obtained for 
the Lynden-Bells' Hamiltonian of type H{q\ 

Incidentally, whenever these microstates of the system in the limit N = 
oo vary sufficiently mildly, then they coincide with the mesostates and map 
bijectively into the macro states. The upshot is: such solutions describe systems 
which in general are truly thermo dynamical] This leads us to tentatively 
propose the following terminology. 

The static maximizers of entropy s(f) over will be called thermostatic 
states. The non-static stationary entropy maximizers will be sanctioned ther- 
mo stationary states; these include some (but not all) of those states which 
are frequently called "rotating equilibrium states." The non-stationary en- 
tropy maximizers are genuinely thermodynamical states. All this applies to 
the N = oo states. 

The finite-A meso- and macrostates for large enough N are then well- 
approximated by these N = oo states and can be categorized in these terms 
as well. This is the content of a law of large numbers as N — > oo over the 
empirical statistics over the particles in an individual system. Note that when 
there is a unique maximizer of entropy, then this law of large numbers and the 
one stated earlier for finite N will give essentially the same results when N is 
large enough, but only then. 



8 Closing remarks 

In this presentation I have tried to explain what equilibrium statistical me- 
chanics is all about. My title "Statistical equilibrium dynamics" summarizes 
this as concisely as I could, in the following sense: Gibbs referred to Boltz- 
mann's stationary ergodic ensembles as being in "statistical equilibrium" (so 
these two words are bound also in my title), but the "dynamics" in my ti- 
tle refers to the meso- or macrostate of the individual members of such an 
ensemble, which as we saw, may very well carry out a nontrivial large scale 
dynamics even in the limit N — > oo. Yet, at any instant of time the individual 
macrostate of the N = oo system is a maximizer of the Boltzmann entropy 
function (-H function) subject to the the conservation laws and other sclero- 
nomic constraints. 

In this presentation I have restricted myself to the general picture, though I 
mentioned specific systems without descending into the details. Those details 
will be given in a joint publication with Carlo Lancellotti |KiLa08j . 
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